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Preface 


This  report  was  prepared  by  the  Coastal  Engineering  Research  Center 
(CERC),  U.S.  Army  Engineer  Waterways  Experiment  Station  (WES),  and  is  a 
product  of  a  reimbursable  study  funded  by  the  U.S.  Army  Engineer  District, 
Los  Angeles  (SPL).  The  investigation  was  conducted  during  the  period 
September  1993  to  January  1994  by  personnel  of  the  Prototype  Measurement 
and  Analysis  Branch  (PMAB),  CERC.  This  report  was  prepared  by  Messrs. 
James  Rosati  III  and  James  P.  McKinney,  PMAB. 

During  the  course  of  the  study,  significant  liaison  was  maintained  between 
WES,  SPL,  the  Ports  of  Los  Angeles  and  Long  Beach,  and  the  Pier  400  design 
team.  Project  management  for  SPL  was  administered  by  Ms.  Jane  Grandon 
under  the  supervision  of  Mr.  Arthur  Shak,  Chief,  Coastal  Engineering  Section. 
Messrs.  Richard  Wittkop  and  John  Foxworthy  were  points  of  contact  for  the 
Port  of  Los  Angeles  and  provided  valuable  input  along  with  their  Pier  400 
design  team  contacts,  Dr.  Kimo  Walker,  Mr.  Russ  Boudreau,  and  Dr.  Paul 
Szwetlot. 

General  supervision  was  provided  by  Dr.  James  R.  Houston,  Director,  and 
Mr.  Charles  C.  Calhoun,  Jr.,  Assistant  Director,  CERC;  direct  supervision  of 
the  project  was  provided  by  Messrs.  Thomas  W.  Richardson,  Chief,  Engi¬ 
neering  Development  Division,  CERC,  and  William  L.  Preslan,  Chief,  PMAB. 

During  the  publication  of  this  report,  Director  of  WES  was  Dr.  Robert  W. 
Whalin.  COL  Bruce  K.  Howard,  EN,  was  Commander  of  WES. 


The  contents  of  this  report  are  not  to  be  used  for  advertising,  publishing, 
or  promotional  purposes .  Citation  of  trade  names  does  not  constitute  an 
official  endorsement  or  approval  of  the  use  of  such  commercial  products. 


Conversion  Factors,  Non-SI 
to  SI  Units  of  Measurement 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
units  as  follows: 


Multiply 

By 

To  Obtain 

feet 

0.3048 

meters 

inches 

2.54 

centimeters 

square  inches 

6.4516 

square  centimeters 

- - — - - * - 

miles  (U.S.  statute) 


1.609347 


kilometers 


1  Introduction 


As  part  of  a  reimbursable  study  for  the  U.S.  Army  Engineer  District, 

Los  Angeles,  a  long  wave  probability  analysis  has  been  performed.  The 
purpose  of  this  work  is  to  supply  the  data  needed  to  predict  the  magnitude  and 
distribution  of  wave  energy  in  the  harbors  after  construction  of  the 
Los  Angeles  Harbor  Pier  400  expansion.  Prototype  wave  data  have  been 
collected  at  offshore  Platform  Edith  as  well  as  selected  inshore  sites  within 
Los  Angeles/Long  Beach  (LA/LB)  Harbors  since  1984  (Figure  l).1  This 
gaging  has  occurred  with  the  support  of  the  Los  Angeles  District  and  the  Ports 
of  Long  Beach  and  Los  Angeles.  Analysis  and  data  collection  methods,  as 


Figure  1 .  Prototype  wave  gage  locations 


I  a  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI  units  is  presented  on 
page  v. 
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i 


well  as  data  availability,  are  documented  in  a  wave  data  summary  report 
(Rosati,  McKinney,  and  Puckette,  in  preparation).  These  wave  data  are  used  to 
document  the  waves  occurring  in  the  harbors  due  to  the  incident  wave  climate. 
The  harbor  response  to  incident  wave  energy  affects  ship  motion  and  therefore 
the  economic  use  of  port  facilities.  This  report  documents  how  the  probability 
analysis  was  performed  and  presents  the  results  for  use  by  the  harbor 
resonance  analysis  performed  by  Seabergh  and  Thomas  (1995). 
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2  Approach 


Transfer  Functions 

The  procedure  for  analyzing  these  wave  data  can  be  viewed  as  a  sequence 
of  transfer  functions.  The  directional  wave  spectrum  outside  the  harbor  is 
transferred  to  short-  and  long-period  waves  inside  the  harbor.  Long-period 
waves  at  each  berthing  site  can  be  used  to  drive  a  ship  motion  model.  Using 
these  transfer  functions  and  long-term  wave  statistics  outside  the  harbor, 
probability  of  exceedance  can  be  computed  for  ship  motions  at  a  berthing  site. 
Prototype  measurements  inside  and  outside  the  harbor  are  used  to  empirically 
estimate  the  wave  energy  transfer  function.  However,  prototype  measurements 
are  available  for  only  a  few  gaging  sites  and  are  representative  of  only  the 
harbor  configuration  existing  while  the  gaging  is  occurring.  Prototype  data  are 
used  to  calibrate  and  verify  the  performance  of  the  model  by  comparing 
prototype  measurements  to  model  measurements  at  the  same  sites  and  for  the 
same  conditions  as  measured  during  prototype  events  (Figure  2). 

The  distorted-scale  harbor  physical  model  is  used  as  the  computational  tool 
to  estimate  wave  energy  at  all  sites  of  interest,  and  for  proposed  new 
configurations.  The  collected  prototype  data  can  only  describe  the  response  of 
the  current  configuration  of  the  harbors  and  therefore  cannot  predict  harbor 
response  when  the  configuration  changes.  Different  expansion  configurations 
were  tested  in  the  LA/LB  physical  model  (Figure  3).  The  existing  configu¬ 
ration  is  referred  to  as  the  base  condition  and  the  proposed  configuration  (the 
only  one  considered  in  this  study)  is  termed  POLA  stage  1. 

Physical  Model  Data 

An  example  of  the  transfer  function  results  obtained  from  the  physical 
model  is  the  amplification  factor  plot  shown  in  Figure  4.  Figure  4  describes 
the  response  of  the  harbor  for  a  given  incident  wave  condition.  In  this  figure, 
if  the  response  is  greater  than  1,  the  waves  inside  the  harbor  are  larger  than  the 
waves  outside  the  harbor  for  the  periods  of  interest.  This  figure  also  shows 
good  agreement  with  prototype  data  measured  at  the  same  location  as  in  the 
physical  model.  The  condition  selected  for  this  study  is  the  uniform  or 
“mean  energy”  spectrum  with  incident  waves  coming  from  the  south.  This  is 
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Long  Period  Wave  Spectra 
Selected  For  Testing 


Figure  2.  Long-period  spectra  events  selected  for  model  testing 

the  most  frequent  wave  condition  incident  to  the  harbors,  but  this  condition  is 
not  necessarily  representative  of  extreme  storm  events.  Physical  model  study 
results  are  described  in  Seabergh  and  Thomas  (1993,  1995).  They  express 
amplification  factors  in  terms  of  wave  height.  This  study  uses  wave  energy 
amplification  factors.  For  this  study,  the  conversion  formula  between  energies 
and  wave  heights  is 


H  = 


(1) 


where  H  is  the  wave  height  and  E  is  wave  energy. 
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Figure  3.  Pier  400  stage  1  model  gage  sites 


Figure  4.  Comparison  of  prototype  and  model  wave  amplification  factors 
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3  Methods 


Platform  Edith  spectral  information  was  used  to  construct  probability 
distributions  of  recorded  wave  data.  The  data  used  for  this  study  covered  the 
time  period  from  1985  to  August  1991.  Although  data  subsequent  to  this  time 
exist,  it  was  most  expedient  for  this  study  to  utilize  these  data  for  which 
previously  computed  intermediate  data  products  exist.  These  data  products 
reduced  the  computational  and  data  retrieval  time  needed  to  perform  this  study. 

Other  workers  (Bowers  (1992),  for  example),  have  defined  the  wave  period 
range  of  30-200  sec  as  the  one  that  primarily  affects  ship  motion  in  harbors. 
The  data  analysis  was  performed  on  two  period  ranges.  One  range  uses  the 
512-sec  and  256-sec  center  periods  from  the  analyzed  wave  spectra,  and 
another  uses  42-  to  170-sec  center  periods.  With  a  spectral  bandwidth  of 
0.00195  Hz,  the  42-  to  170-sec  center  period  range  is  equivalent  to  the  41-  to 
204-sec  period  range  and  the  256-  to  512-sec  center  period  range  is  equivalent 
to  a  204-  to  1,022-sec  period  range.  These  period  range  selections  were  made 
in  order  to  best  describe  the  long  wave  probability  distribution  for  the  harbors 
in  terms  that  are  useful  with  existing  ship-motion  criteria  (Seabergh  and 
Thomas  1995).  Probability  of  exceedance  (PE)  tables  of  wave  energy  were 
constructed  for  both  base  and  stage  1  conditions. 


Data  Processing  Procedure 

a.  Assume  linear  waves  and  transform  each  Platform  Edith  spectrum  by 
multiplying  each  spectral  value  S(f)  by  the  corresponding  amplification 
factor  for  each  model  site. 


He(f)  =  AH(f)S(f) 


(2) 


where 


He(f)  =  harbor  predicted  energy  spectrum 
AH  ( f)  =  model  energy  amplification  factor 
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S  (f)  =  Platform  Edith  energy  spectrum 
/  =  frequency 

b.  Sum  the  energy  belonging  to  each  of  the  period  ranges  of  interest  for 
each  transformed  Platform  Edith  wave  spectrum  ( He ). 


where 

H  S'U  (t)  -  total  harbor  energy  in  period  range 
t  =  center  period  band 

c.  Update  the  tally  of  the  number  of  occurrences  of  each  energy  value  in 
the  bins.  This  is  done  separately  for  each  of  the  two  period  ranges. 
The  number  of  occurrences  Nr  in  a  given  bin  is  equal  to  the  total 
number  of  wave  records  for  which  the  value  of  H  Se  lies  within  the 
given  bin  energy  interval. 

d.  Repeat  steps  1  through  3  until  all  Edith  wave  spectra  have  been 
analyzed. 

e.  Compute  the  probability  of  occurrence  for  each  site,  each  energy  bin, 
and  for  each  period  range. 


po,U) 


*„(0 

N 


(5) 


P02(j) 


Nr2H) 

N 


(6) 
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where 


i  =  index  of  0.35  cm2  energy  bin  associated  with  energy  value; 

j  =  energy  value  associated  with  the  fth  energy  bin  (0.35  cm 

0 

P  0I  2  =  probability  of  occurrence 

N  =  total  number  of  occurrences  (wave  records) 
f.  Compute  the  PE  for  each  harbor  site. 

The  PE  for  energy  was  determined  by: 

PE(J)  -  1  -  £  PO  (Af  -k)  (7 

*=i 


where 


P  O  =  probability  of  occurrence 
PE-  probability  of  exceedance 
e  =  bin  size  (0.35  cm2) 
k  -  bin  index 

g.  Plots  were  constructed  showing  the  base  and  stage  1  PE’s  on  the  same 
page. 

This  method  is  favored  for  several  reasons.  It  minimizes  the  numerical 
roundoff  or  truncation  errors  associated  with  inserting  the  data  into  a  finite 
number  of  probability  of  occurrence  bins  by  performing  this  step  only  once  for 
each  Platform  Edith  spectrum.  It  is  also  the  most  straightforward  realization  of 
the  linear  transfer  function  approach.  A  variety  of  methods  to  compute  the  PE 
distributions  were  explored.  For  example,  multiplying  the  period  range 
average  energy  by  the  averaged  amplification  factors  for  that  range  produces 
very  similar  results. 

A  bin  size  Ae  of  0.35  cm2  was  chosen  in  an  attempt  to  satisfy  the  need  of 
good  probability  resolution  while  representing  the  maximum  energies  with  a 
reasonable  number  (in  this  case  80)  of  bins. 
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Prototype  PE 

Steps  2  through  7  of  the  data  processing  procedure  were  used  on  the  harbor 
wave  gage  data.  These  prototype  PE  plots  are  presented  in  Appendix  A. 
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4  Results 


The  standard  approach  in  both  numerical  and  physical  harbor  wave  models 
is  use  of  a  linear,  frequency-domain  transfer  function  to  transfer  energy  from 
outside  to  inside  of  a  harbor.  Since  energy  in  the  offshore  wave  spectrum  can 
be  low  for  harbor  oscillation  periods,  large  amplification  factors  result. 

Results  presented  in  the  appendices  include: 

a.  For  each  prototype  inshore  site  which  had  a  corresponding  model  site, 
PE  tables  were  constructed  for  the  energy  in  the  two  period  ranges. 

The  PE  results  for  the  prototype  gage  were  plotted  along  with  the  base 
and  stage  1  data  computed  using  the  Platform  Edith  wave  data  and  the 
amplification  factors  from  the  physical  modelling.  These  plots  are 
contained  in  Appendix  A. 

b.  Appendix  B  contains  PE  plots  for  all  model  sites  that  have  both  a  base 
and  stage  1  condition.  These  PE  plots  do  not  include  prototype  wave 
data. 

c.  For  each  model  site,  Appendix  C  presents  a  table  of  the  energy  and 
significant  wave  height  at  the  1 -percent  PE  for  both  the  base  and 
stage  1  conditions.  This  PE  value  was  selected  as  a  comparative 
measure  between  the  base  and  stage  1  conditions.  For  some  model  site 
PE  tables,  80  of  the  0.35-cm2  bins  did  not  reach  the  1-percent  level  of 
wave  energy.  For  these  sites,  the  PE  was  recomputed  using  350  of  the 
0.35-cm2  bins.  Even  this  number  of  bins  was  insufficient  for  some  sites 
to  reach  the  1 -percent  level.  Where  this  occurs,  a  “greater  than”  (>) 
sign  is  placed  next  to  the  maximum  energy  bin  value. 


Chapter  4  Results 


5  Discussion 


It  is  apparent  from  the  plots  and  tabular  results  presented  that  POLA 
stage  1  harbor  energies  tend  to  be  lower  in  most  cases  than  the  base  condition. 
The  plots  in  Appendix  A  show  that  prototype  measured  energies  are  less  than 
predicted  by  the  Platform  Edith  transform  function  approach  (refer  to  Seabergh 
and  Thomas  (1995)  for  a  more  complete  discussion  of  the  harbor  resonance 
analysis).  The  largest  waves  show  the  greatest  disparity  between  the  predicted 
and  measured  PE  distributions.  It  is  important  to  realize  that  the  prototype 
harbor  data  were  measured  under  all  incident  wave  conditions  and  that  the 
model  amplification  factors  used  in  the  predictions  are  for  the  uniform 
spectrum  condition  which  most  represents  the  prevailing  wave  climate.  Large 
storms  tend  to  come  from  a  westerly  direction  that  shelters  the  harbors  from 
the  largest  waves  and  more  directional  spreading  of  wave  energy  occurs  during 
storms  (Seabergh  and  Thomas  1995).  Different  sets  of  amplification  factors 
are  available  from  the  other  long-period  spectra  tested  in  the  physical  model  to 
more  accurately  describe  these  storm  conditions. 

For  engineering  purposes,  the  predicted  energy  values  appear  to  be  conser¬ 
vative:  i.e.,  at  a  given  probability  level,  the  predicted  energies  are  likely  to  be 
greater  than  those  actually  encountered. 


Assumptions 

Seabergh  and  Thomas  (1995)  discuss  some  of  the  assumptions  made  by  this 
analysis  and  why  these  assumptions  appear  to  be  validated  by  the  available 
data: 


...the  long  wave  energy  measured  for  each  frequency  at  the 
prototype  wave  gage  in  the  ocean  is  the  same  energy  transformed 
to  the  harbor  entrance.  ...However  analysis  of  wave  data  indicates 
a  good  correlation  between  wave  energy  at  Platform  Edith  and  the 
harbors.  At  the  harbor  gages  themselves  a  possible  reduction  in 
correlation  of  harbor  response  to  offshore  long  wave  conditions 
could  be  due  to  non-linear  transfer  of  energy  to  other  harmonics, 
but  this  does  not  appear  to  be  a  strong  mechanism  for  a  harbor  of 
large  depth  like  LA/LB  harbors. 
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Additional  Work 


Additional  work  could  be  done  to  reduce  the  statistical  variability 
associated  with  results  and  explore  more  aspects  of  the  collected  prototype 
data.  Additional  work  could  include  using  gages  recently  located  near  the 
entrances  to  the  harbors  as  a  basis  for  determining  amplification  factors.  This 
would  help  to  quantify  any  wave  transformation  effects  between  the  harbors 
and  the  platform.  Classifying  Platform  Edith  and  entrance  gage  data  according 
to  energy  level,  direction,  and  spectral  shape  along  with  other  factors  and 
comparing  these  data  to  the  harbor  wave  data  could  provide  further  insight  into 
the  physical  processes  governing  harbor  response.  Furthermore,  recomputing 
the  results  presented  in  this  report  using  prototype  data  collected  after  August 
1991  would  increase  confidence  in  the  results  due  to  reduced  statistical 
variability.  Higher  order  spectral  analysis  can  help  to  show  the  presence  of 
nonlinear  interactions  in  the  prototype  wave  data. 
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PROBABILITY  OF  EXCEEDANCES 
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FOR  MODEL  SITE  1 


256  -  512  SECONDS 


Appendix  B  Model  Site  PE  Plots 


B3 


PROBABILITY  OF  EXCEEDANCES 
FOR  MODEL  SITE  3 


256  -  512  SECONDS 


PROBABILITY  OF  EXCEEDANCES 
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Table  Cl 

1 -Percent  PE  Values 


Model  Site  Period  Range  Condition  Energy  (sq  cm)  Wave  Height  (cm) 


'0 

Stage  1 

0.593 

3.08 

'0 

Base 

2.529 

6.36 

12 

Stage  1 

1.645 

5.13 

12 

Base 

2.941 

6.86 

Stage  1 
Base 
Stage  1 
Base 

Stage  1 
Base 
Stage  1 
Base 

Stage  1 
Base 
Stage  1 
Base 

Stage  1 
Base 
Stage  1 
Base 

Stage  1 
Base 
Stage  1 
Base 

Stage  1 
Base 
Stage  1 
Base 

Stage  1 
Base 
Stage  1 
Base 


23.640 

42.190 

3.548 

1.368 


34.665 

35.532 

16.107 

49.190 

9.115 

7.473 

23.864 

28.873 


009 

009 

009 

009 

170 

170 

512 

512 

Stage  1 

Base 

Stage  1 

Base 

24.690 

101.340 

1.173 

9.805 

19.88 

40.27 

4.33 

12.53 

010 

170 

Stage  1 

1.785 

5.34 

010 

170 

Base 

10.247 

12.80 

010 

512 

Stage  1 

0.993 

3.99 

010 

512 

Base 

12.664 

14.23 

011 

170 

Stage  1 

0.915 

3.83 

011 

170 

Base 

1.180 

4.34 

011 

512 

Stage  1 

0.646 

3.22 

Oil 

512 

Base 

1.284 

4.53 
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Table  Cl  (Continued) 


1  - 

Model  Site 

Period  Range 

Condition 

Energy  (sq  cm) 

Wave  Height  (cm) 

035 

170 

Stage  1 

5.956 

9.76 

035 

170 

Base 

8.474 

11.64 

035 

512 

Stage  1 

0.348 

2.36 

035 

512 

Base 

1.505 

4.91 

036 

170 

Stage  1 

20.674 

18.19 

036 

170 

Base 

81.915 

36.20 

036 

512 

Stage  1 

3.573 

7.56 

036 

512 

Base 

60.734 

31.17 

038 

038 

038 

038 

170 

170 

512 

512 

Stage  1 

Base 

Stage  1 

Base 

1.329 

1.581 

2.249 

2.615 

4.61 

5.03 

6.00 

6.47 

039 

170 

Stage  1 

3.182 

7.14 

039 

170 

Base 

11.090 

13.32 

039 

512 

Stage  1 

0.348 

2.36 

039 

512 

Base 

5.708 

9.56 

Stage  1 
Base 
Stage  1 
Base 


043 

043 

043 

043 

170 

170 

512 

512 

Stage  1 

Base 

Stage  1 

Base 

66.947 

57.318 

11.557 

6.186 

32.73 

30.28 

13.60 

9.95 

044 

170 

Stage  1 

13.073 

14.46 

044 

170 

Base 

15.579 

15.79 

044 

512 

Stage  1 

4.390 

8.38 

044 

512 

Base 

2.640 

6.50 

045 

170 

Stage  1 

4.798 

8.76 

045 

170 

Base 

6.218 

9.97 

045 

512 

Stage  1 

7.654 

11.07 

045 

512 

Base 

3.504 

7.49 

046 

170 

Stage  1 

9.183 

12.12 

046 

170 

Base 

9.905 

12.59 

046 

512 

Stage  1 

4.440 

8.43 

046 

512 

Base 

2.183 

5.91 
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Table  Cl  (Continued) 


Model  Site  Period  Range  Condition  Energy  (sq  cm)  Wave  Height  (cm) 


Stage  1 
Base 
Stage  1 
Base 

Stage  1 
Base 
Stage  1 
Base 

Stage  1 
Base 
Stage  1 
Base 

Stage  1 
Base 
Stage  1 
Base 

Stage  1 
Base 
Stage  1 
Base 

Stage  1 
Base 
Stage  1 
Base 


75.143 

73.179 

0.696 

0.738 

91.365 

104.315 

7.008 

7.882 

78.582 

86.114 

81.215 

89.265 

24.157 

21.705 

59.515 

59.165 


Stage  1 
Base 
Stage  1 
Base 

Stage  1 


>122.150 

>122.150 

8.332 

10.415 


>44.21 

>44.21 

11.55 

12.91 


057 

170 

Stage  1 

>27.650 

>21.03 

057 

512 

Stage  1 

11.297 

13.44 

058 

170 

Stage  1 

12.904 

14.37 

058 

512 

Stage  1 

9.735 

12.48 

060 

170 

Stage  1 

13.655 

14.78 

060 

512 

Stage  1 

1.880 

5.48 

061 

170 

Stage  1 

5.534 

9.41 

061 

512 

Stage  1 

0.691 

3.33 

063 

170 

Stage  1 

5.973 

9.78 

063 

512 

Stage  1 

0.998 

4.00 
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